Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory.

نویسندگان

  • Ondrej Marsalek
  • Thomas E Markland
چکیده

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects.

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically impro...

متن کامل

Ab-initio study of Electronic, Optical, Dynamic and Thermoelectric properties of CuSbX2 (X=S,Se) compounds

Abstract: In this work we investigate the electronic, optical, dynamic and thermoelectric properties of ternary copper-based Chalcogenides CuSbX2 (X= S, Se) compounds. Calculations are based on density functional theory and the semi-classical Boltzmann theory. Computations have been carried out by using Quantum-Espresso (PWSCF) package and ab-initio pseudo-potential technique. To estimate the e...

متن کامل

Quantum wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems.

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by...

متن کامل

Hydrogen Abstraction Reaction of Hydroxyl Radical with 1,1-Dibromoethane and 1,2-Dibromoethane Studied by Using Semi-Classical Transition State Theory

The hydrogen abstraction reaction by OH radical from CH2BrCH2Br (R1) and CH₃CHBr2 (R2) is investigated theoretically by semi-classical transition state theory. The stationary points for both reactions are located by using ωB97X-D and KMLYP density functional methods along with cc-pVTZ basis. Single-point energy calculations are performed at the QCISD(T) and CCSD(T) levels of theory with differe...

متن کامل

A Theoretical Study on Applying Conformational Analysis of 2-Halo-2-Oxo-1,3,2-Dioxaphosphorinanes

The quantum-chemical calculations on the conformational properties of 2-flouro-2-oxo- 1,3,2-dioxaphosphorinane (1), 2-choloro-2-oxo-1,3,2-dioxaphosphorinane (2) and 2-bromo- 2-oxo-1,3,2-dioxaphosphorinane (3) have been investigated by means of ab initio molecular orbital (HF/6-311+G**) and hybrid density functional theory (B3LYP/6-311+G**) based methods and Natural Bond Orbital (NBO) interpreta...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • The Journal of chemical physics

دوره 144 5  شماره 

صفحات  -

تاریخ انتشار 2016